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Abstract 

We calculate the two-pion light-cone distribution amplitudes in the effective low- 
energy theory based on the instanton vacuum. These generalized distribution ampli- 
tudes describe the soft (non-perturbative) part of the process 7*7 — > vrvr in the region 
where the cm. energy is much smaller than the photon virtuality. They can also be 
used in the analysis of exclusive processes such as 7*p ^ p + 27r, Svr etc. 
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Hadron production in photon-photon colhsions at low invariant masses has become a 
subject of great interest recently [|l|]. Under certain conditions such processes are amenable to 
a partonic description, if the virtuality of the photon, Q^, is much larger than the squared cm. 
energy, W^. The simplest such process, which has intensively been studied 0, is 7*7 7r°, 
which provides a unique framework for studying the pion light-cone distribution amplitude. 
Recently, Diehl et al. have argued that also production of a pair of hadrons, 7*7 hh, 
can be described in factorized form 0]. The amplitudes for these processes contain new 
non-perturbative functions describing the exclusive fragmentation of a quark-antiquark pair 
into two hadronsQ. In addition to the usual parton momentum fraction with respect to the 
total momentum of the hadronic final state, these functions depend also on the distribution 
of longitudinal momentum between the two hadrons, as well as on the invariant mass of the 
produced system, W^. In particular, the authors of Ref.p| consider the production of two 
pions, the soft part of which is contained in a generalized two-pion distribution amplitude. 
In contrast to the one-pion distribution amplitude, which has been studied using QCD sum 
rules and other non-perturbative methods the two-pion distribution amplitude is an 

entirely unknown object. Quantitative estimates of this function are urgently needed for the 
computation of the cross section for such processes. 

In this note we calculate the two-pion distribution amplitude at a low normalization 
point in the effective low-energy theory based on the instanton vacuum 0]. This approach 
has recently been used to study the pion and the photon distribution amplitudes [|], |^. The 
pion distribution amplitude was found to be close to the asymptotic one and consistent with 
the recent CLEO measurements |p. 

By crossing symmetry the process 7*7 hh is related to virtual Compton scattering, 
7*/i '-fh, which can be factorized into a hard photon-parton scattering amplitude and 
an off-forward parton distribution (OFPD), which in addition to the light-cone momentum 
fraction, x, depends also on the longitudinal component of the momentum transfer to the 
hadron, ^. The off-forward isosinglet quark/antiquark distributions in the nucleon have been 
recently been computed in the effective low-energy theory based on the instanton vacuum 
[0. It was found that at x = ±,^/2 the OFPD exhibits characteristic discontinuities^. 
Actually, these would-be discontinuities, which would violate factorization, are reduced to 
sharp but continuous crossovers by the fact that the dynamical quark mass derived from 
the instanton vacuum is momentum dependent. Given the intimate connection between 
the OFPD and the generalized two-particle distribution amplitudes we expect characteristic 
discontinuities also in the latter. We indeed find such would-be discontinuities in the two- 
pion amplitude considered here. However, contrary to virtual Compton scattering, in the 
7*7 — > hh case they do not seem to lead to violations of the factorization assumptions. 

Here we perform a rough model estimate of the two-pion distribution amplitude at a low 
scale (~ 600 MeV), following the approach of Refs.0, ||. Our intention is to discuss qualita- 
tive features of these new objects, such as their dependence on the momentum distribution 

^Similar functions have previously been introduced to describe multi-hadron production at large invariant 
masses [Q. 

^Recently, Radyushkin has discussed the behavior of the off-forward parton distribution near x — ±^/2 
in more general terms using double distributions ||l^ . 
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between the pions, the invariant mass, and isospin, and to calculate the convolution with 
the tree-level hard scattering amplitude. An account of a more extensive investigation will 



be published elsewhere ||TT 



Generalized two pion distribution amplitudes. Our aim is to compute the generalized 
two-pion distribution amphtudes (GDA's), which are defined as 



Hz,C,W') = /dx-e-t^^^^"(7r>i)vr''(p2)l^(a:)nT^(0)|0) . (1) 



4tt 



=0, x±=0 



Here, n is a light-like vector (n^ = 0), which we take as = (1,0,0, 1). For any vector, 
V, the "plus" light-cone coordinate is defined as = {n ■ V) = + V^; the "minus" 
component as V~ = — V^. The outgoing pions have momenta Pi,P2, and P = pi + p2 
is the total momentum of the final state. Finally, T is a fiavor matrix (T = 1 for the 
isosinglet, T = for the isovector GDA). For brevity we suppress the isospin indices a, h on 
the function $. 

The generalized distribution amplitude, Eq. (Q), depends on the following kinematical 
variables: the quark momentum fraction with respect to the total momentum of the two-pion 
state, z\ the variable C, = pf/P^ characterizing the distribution of longitudinal momentum 
between the two pions, and the invariant mass (cm. energy), W"^ = P^. 

In what follows we shall work in the reference frame where Pi = 0. In this frame 

P- = PI = pI - = = w\{i-0-K. (2) 

From these relations one obtains the following kinematical constraint: 

2 

C(l-C)>:^. (3) 

We shall compute the GDA, Eq. d), using the effective theory of pions interacting with 
massive "constituent" quarks, which has been derived from the instanton model of the QCD 
vacuum 0. The coupling of the pion field to the quarks, whose form is restricted by chiral 
invariance, is described by the action 



Sint = / d^x %l){x)^M{d^) U^^{x) ^M{d^)tlj{x), (4) 



where 

U^^{x) = e'^«""'^"(")/^- = 1 + — 757r"(a:)r<^ - AT7r'^(a;)7r'^(a;) + .... (5) 

Here, F^^ = 93 MeV is the weak pion decay constant 

The momentum-dependent dynamical quark mass, M{p), plays the role of an UV regula- 
tor. Its form for Euclidean momenta was derived in Ref. 0. The momentum dependent mass 
cuts loop integrals at momenta of order of the inverse average instanton size, ~ 600 MeV. 
One should note that the value of the mass at zero momentum, M(0), is parametrically small 
compared to p~^; the product M(0)p is proportional to the packing fraction of the instantons 
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in the medium, which is a small parameter fundamental to this picture. Numerically, a value 
M(0) = 345 MeV was obtained in Ref.0. 

Computation of GDA in the effective chiral theory. With the simple action Eq. @ we can 
proceed to compute the matrix element defining the GDA, Eq. (|I]). In the leading order of 
the 1/iVc-expansion the result for $(2;, (, W^) is given by the sum of the two Feynman graphs 
shown in Fig. |I]. [The quark propagator here contains the dynamical quark mass.] The first 
graph. Fig. |l|a, evidently contributes only to the isosinglet part of the GDA (T = 1), whereas 
the second graph. Fig. ^3, contributes to both the isosinglet and isovector part (T = r^). 

We first consider the contribution from Fig. |l|a. By straightforward calculation we obtain 

A-i N P+ r rl'^h I 

<|.«(z,C,l^^) = ^^^5^'Tt{T) j ^^5{k^ - zP^)^ M{k)M{k - P) 

zM{k-P)-{l-z)M{k) 



[P - M{ky + iO] [{k - Pf - M{k - Py + zO] 

From this we immediately see that ^^^\z, W"^) is independent of C, and also that 

= z,W^), (7) 

a property which actually follows from C-invariance [Q. 

Let us first evaluate the integral in Eq. (^) neglecting for a moment the momentum 
dependence of the constituent quark mass. In this case the result takes a simple form, 

^m(zw^) = _ ^^^o tt(t) eizd- z)](2z -i) f ^^^^ - 

(8) 

where Mq = M(0) and 6 denotes the step function. This expression is nonzero at the 
endpoints, z ^ and 1. Such behavior of the distribution amplitude would violate the 
factorization theorem for the process 7*7 tttt, because the amplitude for this process 
contains the integral 

new') = £dz^^^^z,CW'), (9) 

which would be divergent if $ were nonzero at the endpoints. However, such conclusion 
would be premature, as one can easily see that the momentum dependence of the constituent 
quark mass becomes crucial at the end point (see the discussion of this point in Refs.[|T^, ^). 
This means that when computing the isosinglet GDA one cannot neglect the momentum 
dependence of the constituent quark mass. For the numerical estimates we can employ a 
simple numerical fit to the momentum dependent constituent quark mass obtained from the 
instanton model of the QCD vacuum. 
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Results of a computation of the distribution are shown in Fig. ^ where we have plotted 
the contribution ^^^^ to the isosinglet GDA as a function of z at various values of W"^. 
[Remember that this contribution is ^independent] . It is useful to expand the isosinglet 
GDA in Gegenbauer polynomials of index 3/2, which are the eigenfunctions of the evolution 



kernel 13 



odd 



Due to the antisymmetry with respect to z ^ 1 — z only polynomials of odd degree appear 
in the expansion of the isosinglet part. The numerical results for the first few coefficients 
can in the momentum range < < 4Mq be approximated by the formsQ 



0.5556 



0.75 
0.036(1 - 6.2 



0.82 



,^^,2^ 0.0036 
a^iW^) = —Tj^, (12) 




where W"^ is taken in GeV^. One sees that the coefficients of the expansion Eq. (|TTD decrease 
rapidly with increasing the order of Gegenbauer polynomials. We also give the result for the 
contribution of to the isosinglet part of the integral of Eq. (^). With good accuracy the 
numerical results can be fitted by 

/«(iy2) = -rtr(T)^^. (13) 

^ ~ To 

Let us now turn to the second contribution to the GDA, given by the graph in Fig. 
One can easily convince oneself that the neglection of the momentum dependence of the 
constituent quark mass in this case does not lead to any violation of factorization. There- 
fore, in order to make the discussion more transparent we shall in the following neglect the 
momentum dependence of the constituent quark mass in the contribution of graph Fig. [l|b. 
The logarithmic divergence of the loop integrals can always be absorbed in the pion decay 
constant, F.^, which is defined by the logarithmically divergent integral (written here as an 
integral over Euclidean momenta) 
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= AN, I ^ (14) 
' ' (27r)4(A;2 + M2)2" ^ ^ 



■^The functions given here should be considered as purely numerical fits. In particular, the powers in the 
denominators have no physical meaning, but are simply chosen such as to represent the numerical results in 
the most compact form. 
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All formulae below can easily be generalized to the case of a momentum-dependent mass. 

Computing the Feynman integral corresponding to graph Fig. |l]b we obtain for this 
contribution to the GDA 
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where 



+Tr(T) 5' 



.<C M^N, f (Pk 



ab 



(15) 




(27r)2 kl + Mi - W^z{l - z) 

k±-p± + m\z - W^z{l - C) 



(16) 



Mi)C - 2k^ ■ p^z - mlz + W^{1- C)z^ 
the value at 2; > C can be obtained from the above expression as 

<P^'\z,CW') -0(2)'(^<)(l-z,l-C,l^^). (17) 

Some comments are in order here. The expressions for the GDA given by Eqs. (§) and 
( |l5i |16|, |l^ are valid in the parametrically wide region of invariant cm. energies 4m^ < 



W"^ < 4Mq . For W"^ > 4Mq the graphs of Fig. |l| develop an imaginary part which is related 
to the presence of a two-quark threshold in this model. This singularity is not physical, 
because for momenta of order ~ 600 MeV degrees of freedom not accounted for in the 
effective theory {e.g. heavy resonances) start to play significant role. These contributions 
can in principle be estimated phenomenologically by adding the coupling of the resonances 
to constituent quarks to the effective action, Eq. @. Note also that the GDA's computed in 
this model are real at low only in the leading order of the 1/A^c^expansion; an imaginary 
part appears in the next-to-leading order due to pion loop contributions, which can be 
estimated using chiral perturbation theory techniques. 

The contribution to the GDA given by Eqs. (^ |1^, |1^ exhibits discontinuities ai z = ( 
and z = 1 — (. The nature of these discontinuities is the same as those found in the non- 
diagonal parton distributions computed in this model |12[. As was shown in Ref. |[L2||, the 



inclusion of the momentum dependence of the constituent quark mass smears this discontinu- 
ity. However, in contrast to the case of non-diagonal parton distributions, the discontinuities 
in the two-pion GDA do not lead to violation of factorization, so we can be negligent and 
evaluate Eqs. (^, |I^) with a constant quark mass. 

Let us study the result Eq. (|^) in the chiral limit, m,r = 0, and expand in powers of W^. 
In this way we obtain an analytic expression for the GDA, 

,-HM,w^) .{1 + ^^iMi^ip^ + ...}, 

,^^.(.C, W^) -(1 - g^ii^ilzOi^ ....). ,18, 
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From this expression we find for tlie lowest two moments 



ldz¥'\z,CW') = 2(2C-l)Tr(T[r^r^])|l + ^^^y^ + ...|, (19) 
jyz{2z-l)^'^'\z,CW') = rtr(T)|i-2C(l-C) 

Tlie first moment, Eq. ([T9|), corresponds to the expansion of the pion electromagnetic form 
factor for small time-like momenta. The pion charge radius read off from Eq. (|19D is 



{r'%.^. = (21) 

which coincides with the result obtained previously in Ref. 0. Note that only the isovector 
part of the GDA contributes to the first moments, as it should be on grounds of C-parity. 

Of $(2) only the isosinglet part contributes to the integral Eq. (^. This contribution to 
the integral is given by 

I^'\C,W') = 4r^Tr(T) |-(l + ilog[C(l-C)]) 

+ ^^(4C(l-C) + 3C^log[C(l-C)] + 3(l-2C)log(l-C))}. (22) 

Note that the R.H.S. is symmetric with respect to C — > 1 — C- 

The total result for the integral Eq. {^) is given by the sum of Eqs. (p!3D and (^)). In 
Fig. ^ we plot the total result as a function of ( at several values of W^. We see that the 
absolute value of this integral for intermediate values of ( increases with W'^, and that the 
dependence on ( is strong only for small values of ({1 — (), a region which is difficult to 
access because of the kinematical constraint, Eq. (^). 

Finally, in Fig. ^ we show the total result for the isosinglet GDA, given by the sum of 
Eq. (D, and Eq. dH), for a value of C = 0.25. Note the different symmetry of the 
contributions with respect to z 1 ~ z, as can be seen also from the analytic expressions. 
We emphasize once more that the discontinuities in ^^'^^ are smeared out when taking into 
account the momentum dependence of the dynamical quark mass. 

Asymptotic behavior of the GDA. The scale dependence of the GDA is governed by the 
usual Efremov-Radyushkin-Brodsky-Lepage evolution |TE|, |^. This is obvious when one 
notes that the hard-scattering kernel of the amplitude for 7*7 tttt is the same as that for 
a process with only one meson in the final state (with appropriate quantum numbers). The 
evolution equations can thus be solved, in analogy to those for the usual pion wave function 
[0], by expanding in eigenfunctions of the evolution kernel. In the isosinglet part of the 



GDA one has to take into account the mixing with gluon operators the evolution of the 
isosinglet GDA calculated in the effective low-energy theory will be considered elsewhere 
[ p!l|] . However, one can easily establish the scale dependence of the isovector part of the 
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GDA, which does not mix with gluons. This function does not enter in the amphtude for 
7*7 — i> vrvr, but it can appear in exclusive processes of the type 7*p — ^ Xp, where X is a 
two- (or more) pion state, cf. the discussion below. The asymptotic behavior of the isovector 
GDA is given by 

HzX.W) = 3 Tr (T[r", r"]) ^(l-^)(2C-l)F.(W^^), (23) 

where FT^iW"^) is the pion electromagnetic form factor in the time-like domain. 

Conclusions. In this note we have given a first model estimate of the two-pion GDA's at a 
low scale. A number of qualitative features of the two-particle distribution amplitudes have 
emerged, which are consequences of the general structure of the effective low-energy theory 
and should thus be general, even given the crudeness of the quantitative approximations 
made. 

First, the GDA's generally exhibit discontinuities at z = ( and 1— C, which are in complete 
analogy to the discontinuities found in the off-forward parton distribution functions^ |T2| . 
We have also noted that, as in the OFPD's, the momentum dependence of the dynamical 
quark mass derived from the instanton vacuum turns these discontinuities into continuous 
but sharp crossovers. However, since the discontinuities of the GDA's are not at odds with 
the factorization of the amplitude for 7*7 — >■ hh we have been negligent about implementing 
the momentum-dependent mass and estimated the amplitude with the discontinuous GDA's. 
Second, we have found that the isosinglet and -nonsinglet GDA's for the two-pion final state 
computed at the low scale show different symmetry properties with respect to z ^ 1 — z. 

Our results concerning the VT^-dependence of the GDA's and the total amplitude should 
be regarded as crude estimates. In particular, the imaginary part of the functions at W"^ > 
4Mq should not be regarded as physical, since it occurs in a region where the effective theory 
is no longer applicable. A more conservative approach would be to calculate the GDA's in 
an expansion in the external pion momenta and the pion mass, adopting, so to speak, the 
philosophy of the chiral Lagrangian. The physical imaginary part, which is due to final state 
interactions, appears in the next order of the 1/Nc expansion and can be estimated using 
the methods of chiral perturbation theory. This approach will be developed in a forthcoming 



paper |Tl| 



We note that the studies reported in this note can easily be generalized to distribution 
amplitudes for n > 2 pions in the final state. Such multi-pion GDA's play a role in the anal- 
ysis of exclusive processes of the type 7*p — > Xp, where X = 27r, Stc ||15[. For instance, 
the two-pion isovector GDA computed here can be used to analyze the 27r background in 
the diffractive leptoproduction of rho mesons. 
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''in Ref.|iJ] only the OFPD's of the nucleon were computed. The OFPD's in the pion at the low scale 
exhibit would-be discontinuities similar to those in the nucleon. A complete study of the OFPD's in the 
pion will be published shortly jlll . 

^We are grateful to L. Frankfurt and M. Strikman for discussion of this point. 
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(a) (b) 

Figure 1: The diagrams defining the two contributions to the two-pion distribution am- 
phtude in the effective low-energy theory, and The solid line denotes the quark 

propagator with the dynamical quark mass, the solid points the vertices defined by Eq. (^). 
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Figure 2: The contribution (^^^\z,W'^) to the isosinglet two-pion distribution amphtude, 
Eq. (|^), as a function of z for various values of W^. This contribution is independent of (. 
Solid Une: W = 0.3 GeV; dashed Ime: W = 0.5 GeV, dotted Une: W = 2Mo = 0.69 GeV. 
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Figure 3: The integral determining the 7*7 —>■ tttt amphtude, Eq. (P), given by the sum of 
Eqs. ([T3| ) and (p2D, as a function of ( for various values of W^. Solid line: W = 0.3 GeV; 
dashed line: W = 0.5 GeV, dotted Une: W = 2Mo = 0.69 GeV. 
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Figure 4: The total result for the isosinglet two-pion distribution amplitude, ^{z, (, W^), as 
a function of for ( = 0.25 and = 0.25 GeV^. Dashed line: contribution Eq. (j^); 
dotted line: contribution Eqs. ([15| , |T6|, solid line: total. 
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